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Abstract 

Motivated by the recent contradiction of the superconducting pairing symmetry in the angle- 
resolved photoemission spectra (ARPES) and the nuclear magnetic resonance (NMR) data in the 
FeAs superconductors, we present the theoretical results on the phase diagram, the temperature de- 
pendent Fermi surfaces in normal state, the ARPES character of quasiparticles and the spin-lattice 
relaxation 1/Ti of the two-orbital t-t -J-J models. Our results show that most of the properties 
observed in iron-based superconductors could be comprehensively understood in the present sce- 
nario qualitatively, indicating that the pairing symmetry of the ironpnictides is anisotropic nodeless 
d x 2_ Jiy 2+S x 2 y 2-wave, mainly originating from the band structures and the Fermi surface topology. 

PACS numbers: 74.20.Rp, 71.7Q.Ej, 74.25.Dw 
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I. INTRODUCTION 



Following the discovery of the first high temperature superconductivity in the copper- 
based compounds two decades ago a second class of high temperature superconductors 
has been recently reported in iron-based pnictides 0, 3, 4, 5, 6, 7|, in which the transition 



temperature T r ca n be as high as 55K 
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251 ] efforts have been devoted to 



understand the nature of new superconductors. Resemble to cuprates, layered iron pnictides 
consist of the conducting FeAs layers and the ReO layers (Re represents the rare earth 
elements, such as La, Ce, Pr, Nd, Sm, and etc) which provide carriers to the FeAs layers. 
The antiferromagnetic spin density wave (SDW) order is suppressed upon substituting a 
few percent O with F, and the new compounds become the superconductor (SC) below 
T c . On the other hand, there are a few of considerable differences between cuprates and 
ironpnictides. For example, the first difference is that the undoped iron-oxypnictides with 
the SDW order exhibit poor metallic conduction, rather than the Neel AFM insulator in 
undoped cuprates; the second one is that the spin and magnetic moment of Fe are much 
smaller than expected, completely different from the Cu spin in cuprates, etc. These facts 
suggest that there exists significant difference between the ground states of ironpnictides 
and cuprates. 

First-principle electronic structure calculations provide the first evidence of the difference. 



It has been shown 



23, 



29 



30 



3l| that in LaFeAsO, the most of the spectral weight 



close to the Fermi energy is contributed from the Fe 3d orbital, and the Fermi surfaces (FS) of 
LaFeAsO consists of two hole-type circles around T point and two electron-type co-centered 
ellipses around M point [32j. This implies that the multiorbital character is dominant in the 
undoped FeAs superconductor, contrast to the single-orbital character in cuprates. Up to 
date, many tight-binding multi-orbital models have been proposed to reproduce the Fermi 
surface (FS) character and the band structures near Fermi level Ep in ReFeAsOi_ x F x and 
Bai .T-K.T-Fe^As-?. To this end, once the electron-phonon mechanism is precluded by Boeri et 



al. 



331 ] . the electronic mechanisms are proposed as the driving force of the SC paring, such 



as the Coulomb interaction between Fe 3d electrons [22j, [23|, [2j3|, |34J or the antiferromagnetic 
exchange coupling between the nearest neighbor Fe sites and the next nearest neighbor site, 

33], three-orbital ones [23], and four-orbital ones 



such as the two-orbital models 



35 



36 
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381 ]. even the five-orbital model 2l|, |39[]. Among these models, the two-orbital t-t'-J-T 



model 



35 



361 ] is minimal, and can reproduce the key characters of ironpnictides, such as the 



complex Fermi surface and the multi-orbital degeneracy of the Fe 3d electrons [3a, [36] and 
the band structures near Ep, as well as the stripe antiferromagnetic or spin- density- wave 
ground state. 

In this newly discovered superconductor, the most intriguing issues are the pairing mech- 
anism and pairing symmetry. To date, various possibilities of the SC mechanism and pai ring 



symmetry for Fe— pnictide 
and experiemntally j^ll . 



(i-wave 
wave 



/e 3 



superconductor are proposed theoretically [22], 



23 



25 



34, 



35 



42 



22 



40] 



43L Q, 3- The pairing symmetries range from spin singlet 



47 ] to s-wave 25|, [34], or the mixture of S x 2 y 2 and d x 2_ y 2 [40], or spin triplet p- 



231 ] . At present, the nature of the SC gap observed experimentally is very different 



rom authors to authors and from samples to samples, ranging from one gap 9| to two gap s 



43 



44" 



481 ] . and from isotropic [42j, |43j or anisotropic fully gap [41] to line node gap [44]. 



On the one hand, the angle-resolved photoemission spectroscopy (ARPES) measurements 



41 



42 



431 ] . a direct measurement to the quasiparticle spectra and the SC gap, observed 



that all the gaps are nodeless around their respective Fermi surface sheets, which is different 
from the situation of the cuprates. On the other hand, however, the spin-lattice relaxation 
rate in the nuclear mag netic resonance (NMR) experiments strongly suggest the existence 



of nodes in the gap 



44j . which share partial features of cuprates and MgB 2 . The question 



thus arises whether the contradiction between the ARPES and the NMR data is due to the 
sample qualities or it reflects the two sides of the new superconductors ? 

In this paper, we start from the minimal two-orbital t-t'-J-j' model j40j, which has the 
same topology as the band structure of the iron based superconductors, and obtain the 
mean-field phase diagram of the extended t-t'-J-j' model, the quasiparticle spectra and its 
ARPES manifestation, and the spin-lattice relaxation rate 1/Ti. Our results demonstrate 
that the various pairing symmetries are stable in the mean-field phase diagram, and the 
T 3 -like behavior in the spin-lattice relaxation rate may coexist with the anisotropic nodeless 
SC gap, as observed in ARPES and NMR experiments. 
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II. MODEL HAMILTON! AN AND METHOD 



Based on the band structures results and theoretical analysis, the twofold-degenerate 
dxz/dy Z orbits are essential for the ironpnictide superconductors. We depict such physical 
processes with the two-orbital Hubbard model, 



<ij>a/3a <^ij^>a/3a 

ia iaa' 

-J H cllAlv'Cva'CMv + J H ^ Ct^'tC^ia' | (1) 
icrcr' iaj^a' 

here c\ aa creates a d xz (a=l) or d yz (a=2) electron with orbital a and spin a at site Rj. 
t and t' denotes the hopping integrals of the nearest-neighbor (NN) and the next-nearest- 
neighbor (NNN) sites, respectively. U, U' and J# are the intra-orbital, inter-orbital Coulomb 
interactions and the Hund's coupling. Formally, the t-t'-J-j' model can be derived from 



Eq.(l) in the atomic limit 49j, although far from strict. The t-t'-J-j' model reads, 

H = H t _ t , + Hj_y (2) 

with the kinetic energy term, 

Ht-t' = ^y^tkxz — / i ) c ll - C fclff + ( £ kyz ~ ^) c k2a C k2a + ^kxy{4\a C k2a + C k2a C kla)] 
ka 

with the notations 

£kxz — — 2ti cos k x — 2t 2 cos k y — 4t 3 cos k x cos k y , 
£kyz — — 2^2 cos k x — 2ti cos k y — 4t 3 cos k x cos k y , 
£kxy = -4t 4 sinfc x sin k y . 

Here the components of t*f are t^, 1 =t 1 =-l and t x 2 =t 2 = 1.3, and these of t'" 13 are t 3 = t 4 = 



—0.85 |40j. The electron is filled with n=2 (or half-filling) on a two-dimensional square 
lattice. The interaction term reads, 

H j-f = J (^ ia ' ^ ~~ 1 nia ' n iP) + J ^2 (^ia ■ Sj/3 - -n ia ■ Ujp) (3) 

<ij>a/3 «ij»a/3 

where J and j' are the NN and the NNN antiferromagnetic couplings, through the concrete 



values of J and j' are still in debate. Haule et al. estimated that J / j' ~ 2 



In this 



/ — * 

paper, we assume J=0.7 and J =0.3. Sf is the spin operator of the electron in the a-orbit, 
and riia is the particle number operator. The orbital indices a and f3 run over 1 and 2. 
Introducing the following order parameters, 



A 1 ? = 

x/y 


J < C lm\ C ]ml >' U 


= i±x/y) 




A 2m 

x±y 


J < C lm1 C ]mi >> U 


= i±(x±y)) 




pm 

x/y 


^ C imtr C jmtj >) (j ' — 


i±x/y) 




pm _ 
^3 — 


^ C ima C jma >j (j = 


i ± (x±y)) 


(4) 



the interaction term in Eq.(4) could be decoupled within the framework of the mean- field 



approximation 



401 ] . Slightly different from Seo etal. mean-field ansatz, we consider the 
contributions of the kinetic order parameters P x / y and P3, since the itinerant character of 
the Fe 3d electrons is considerable. So one could obtain the mean-field Hamiltonian and the 
ground state energy. Minimizing to the ground state energy gives rise to the self-consistent 
equations for the order parameters A and P, 

k 

J' 

A *±V = nYI C ° S ( k * ± k v"> < C -kmAml > 
k 



2N 

k 



P™ = ^^2cos{k x )cos(k y ) < n kma > (5) 

ka 

The inter-orbital pairing parameter < cj^cj^ > is neglected due to its very small value 
[4^]. This implies that the p-wave pairing symmetry is precluded. From Eq.(4), one notices 
that for each orbital, the hopping along the x-direction is not equivalent to that along the 
y-direction. This leads to, for example, A^. m ^A^ m . Due to the equivalence between the d xz 
and the d yz orbits, the order parameters of the two orbits are symmetric once interchanging 
the x-direction and the ^/-direction. Therefore, throughout this paper, we mainly focus the 
order parameters in the first orbit. The results on the second orbit can be arrived at if 



rotating the to the y-axis [401 ] . 



Within the present scenario, we could obtain not only the groundstate phase diagram, but 
also the quasiparticle spectra in the normal and the SC states. The temperature dependence 
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FIG. 1: (Color online) The phase diagram of the t-t'-J-j' model at doping concentration x=0.18 
from the half-filling. N denotes the normal state, the other four phases are SC 

of the Fermi surface in normal state and that of the spin-lattice relaxation rate in the SC 
state can also be obtained. Among these quantities, the spin-lattice relaxation rate in the 



NMR experiment is expressed as 
Tin 2 



J J N s (E)N s (E')f(E)[l - f(E')6(E - E')dEdE'} 



(6) 



T u k B T 

Providing 1/Tin in the normal state satisfies the Korringa law, the spin lattice relaxation 
rate 1/T\ a becomes jsij: \/T\ a oc (ksT) ■ T\n/Tx s . 



III. THEORETICAL RESULTS AND DISCUSSIONS 
A. Phase diagram 

Different from Seo et aZ.'s phase diagram, we find five stable phases in our phase diagram, 
which are shown in Fig.l. The first one is the normal phase, denoted N in Fig.l. Obviously, 

6 



when the superexchange coupling J and J' are too small to provide the SC pairing glue, 
the kinetic energy is dominant, and the electrons stay in the normal state. Among the 
SC phases mediated through the spin exchange couplings, large J and small J' favor the 
S X2+m2 (oc cos(kx) + r]cos(ky))(r] > 0) SC symmetry, which is the combination of the S Xp+ y 



SC symmetry and d X2 _ y2 SC symmetry. Here r\ is a constant. Different from Seo et al. 40j . 
it is not a mixed state of the S X2+y2 SC state and the d X2 _ y2 one, since the relative phase of 
the coefficients of the S X2+y2 component and the d X2 _ y2 one is fixed, hence it is a pure state. 
On the other hand, small J and large J' favors the S X2V2 oc cos(kx + ky) + cos(kx — ky) 
SC phase, since the NNN hopping integrals along the x-axis and the y-axis are identical. 
As seen in Fig. la, in the regions where J and J' compete with each other, the third 
SC phase appears as the combination of the S X2+vy2 and the S X2V2 symmetries. For the 
case in this region with J = 3 and J = 1.5, the SC order parameters are A^ 1 =0.241, 
A^ 1 =0.555 and A^™^=0.060. And the fourth SC phase appears with the combination of the 
d X2 - m2 ((x cos(kx) —rjcos{ky)){rj > 0) and S X2V2 symmetries. For the case in this region with 
J=0.7 and j'=0.3, the SC order parameters are A^-0.039, A^=0.029, and A^=0.012. 
The components of the SC order parameters strongly depend on the relative magnitude of 
j' with respect to J in different SC regions. As we can see the ratio A^/A* 1 =R in Fig. lb, 
the change of the ratio R is discontinuous with the increase of j' at fixed J, implying that 
the quantum phase transitions between these SC phases are the first order. 

B. ARPES and T-Dependence of SC Gaps 

The angle-resolved photoemission spectra (ARPES) experiment provides direct informa- 
tion about the quasiparticle spectra in normal state and the symmetry of the SC gaps in 
the SC phase. We present the SC gap character of the t-t'-J-j' models for the parameters in 
LaFeAsO in Fig. 2. Notice that the benchmark of the Fermi surface adopted here to plot the 
SC gap is at T = 0.2. Our results show that the SC gap structure exhibits fourfold symme- 
try, as seen in Fig. 2a, in accordance with the Fermi surface topology of the present model. 
The amplitudes of the SC gaps on the different Fermi surfaces are different, suggesting the 
nature of a two-gap SC. The SC gap structures are anisotropic and nodeless, as shown in 
Fig. 2a. Though the SC gap structure likes an anisotropic s— wave pairing symmetry, it is in 
fact the symmetry of the d T 2_ m ,2 +S x 2,,2-waye. Our result reaches an agreement with recent 



ARPES experiments by Kondo etal. 4l|. We also notice that the relative variation of the 
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FIG. 2: (Color online)The angle dependence of the SC gaps near the small hole FS (red line) and 
the large hole FS (orange line) around the T point in the polar coordinates. Theoretical parameters: 
J = 0.7, j' = 0.3, and doping concentration x=0.18. 



anisotropic gaps 



is less than 25 %, so our results do not conflict with Ding et al.'s 43| and 



Zhou et a/.'s [42j reports. 

Fig.2b shows the temperature dependence of the SC energy gap A along the 9 = 
and 9 = 90 deg. With the increasing of the temperature, the SC order parameter 
in the small hole-like FS sheets decrease monotonously, as observed in the ARPES 
experiments 



42 



431 ] . Near T = 0.04, both the SC gaps exhibit a small dip, which 



is attributed to the interactions among the three SC energy scales. Obviously, the 
nodeless and the anisotropy of the SC gaps originate from the unique gap structure, 
A m = A^j? y (cos k x cos k y ) + ( Aj. m cos k x — A^ m cos/c y ). And such a symmetry is also in 
agreement with the analysis on the SC pairing symmetry 52j, |53|]. However, one finds 
that the magnitude of the gap in the large FS is bigger than that in small FS in the T 
point, opposite to the ARPES result. This may be due to the tight-binding parameters of 
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FIG. 3: (Color online) (a) The Fermi surface for the different temperature, T=0.17, 0.2, 0.5. (b) 
The theoretical parameters are the same as the Fig. 2. 

the two-orbital t-t'-J-j' model, which only describes the topology structure of the Fermi 
surfaces of the FeAs superconductors, and does not contain all the details of the Fermi 
surfaces and the band structures in ironpnictide compounds. Hence, we expect that the 
refined tight-binding parameters will improve our results in the further study. 



C. T-dependence of Fermi Surface 

The temperature evolution of the FS in the normal state is shown in Fig. 3, when T 
is higher than the critical temperature of the SC phase transition. The two hole-like and 
electron-like FS sheets can be clearly identified around the T and the M points at x= 
0.18, as observed in the ARPES experiments and the first-principle electronic structures 
calculations. Interestingly, the hole-like FS sheets expand a little with the increasing of 
the temperature. In contrast, the electron-like FS sheets shrink very acutely. This may 
demonstrate that the electron-like FS sheets are more important for the occurrence of the 
SC state, in consistent with Dai etal. prediction. 221 ] . 



D. Spin-Lattice Relaxation Rate in NMR 

Although many experimental measurements, such as the Andreev reflection 



exponential temperature dependence of the penetration depths [55( and the ARPES [42j, |43] 



5JL the 
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FIG. 4: (Color online) Temperature dependence of the spin lattice relaxation rate in the t-t -J-J 
model. The red arrow indicates the SC critical temperature T c , The red line is the T 3 law for 
comparison. Inset shows the detail near T c . The theoretical parameters are the same as in Fig. 2. 

observe the nodeless gap function in the SC phase of ReFeAsOi_ a; F :c and Bai _. r K r Feo As? 
compounds, the line nodes in the SC gap was also suggested by the NMR experiment 44]. 
The two characters in the NMR experiment supported the line nodes: lack of the coherence 
peak and the T 3 behavior in the nuclear spin-lattice relaxation rate, 1/Ti. Using the gap 
function obtained in this paper, we calculate the spin lattice relaxation rate 1/T S , and 
the numerical result is shown in Fig. 4. We also plot the T 3 law (the red line) for a 
comparison. It is found that over a wide temperature range, the spin lattice relaxation rate 
in the present model can be fitted by the T 3 law, in agreement with the observation of the 



NMR experiments 



44 



56(]. 



A small coherence peak appears around the critical transition temperature, as clearly 
seen in the inset of Fig. 4. Experimentally, such a small coherence peak may be easily 
suppressed by the impurity effect or the antiferromagnetic spin fluctuations, similar to the 
situations in cuprates. This leads to the missing of the Hebel-Slichter coherence peak in the 
NMR experiment in ironpnictide SC. With the decreasing of the temperature, one find a 
drop in the s pin lattice relaxation rate, 1/T ls , consistent with the observation of the NMR 
experiments [44j. Such a behavior deviating from the T 3 law may contribute from the multi- 
gap character of this system, and such a drop reflects the different SC gaps in different orbits. 
Surely, more meticulous studies are needed in near future. We also notice that Parker etal. 
bund the extended s± SC gaps also can give the same NMR relaxation rate in SC pnictides 



571 ] . The inter-band contribution to the spin-lattice relaxation deviating from the T 3 law 
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was suggested by Parish et al. [58 



Therefore, we find that the numerous unusual properties in the normal state and the SC 
phase of newly discovered FeAs superconductors could be interpreted in the t-t -J-J model, 
showing that this model is a good approximate model to describe the iron-based supercon- 
ductors. Within this scenario, the superimposition and mixing of the pairing electrons in 
the two orbits contributes to the anisotropic nodeless SC pairing symmetry. Such a pairing 
symmetry assembles the characters of usual (i-wave and s-wave, hence shares the properties 
of the usual <i-wave SC, like cuprates, and the s-wave SC, such as MgB 2 ^2]. Nevertheless, 
to compare the theoretical results with the experimental observation quantitatively, more 
subtle band structures in the t-t -J-J model are needed, though the present two-orbital 
model gives the qualitatively correct properties. Also the present constrained mean-field 
approximation should be improved in further study. 



IV. SUMMARY 

In summary, starting with the t-t'-J-J model, we obtain the mean-field parametric phase 
diagram at the doping concentration x=0.18, and find the normal state and a new SC phase 
in the phase diagram, different from the literature. With the decrease of the temperature, a 
nodeless and anisotropic d x 2_ rjy 2+S x 2 y 2-wave gap structure emerges. But due to the multi- 
gap character, the anisotropic gaps open on the hole Fermi surfaces and the electron Fermi 
surfaces are different. The T 3 law of the spin lattice relaxation rate, 1/Ti, is also interpreted 
in the present t-t -J-J model. 
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